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Abstract. We consider the motion of a non relativistic quantum particle in 
subject to n point interactions which are moving on given smooth trajec- 
tories. Due to the singular character of the time-dependent interaction, the 
corresponding Schrodinger equation does not have solutions in a strong sense 
and, moreover, standard perturbation techniques cannot be used. Here we 
prove that, for smooth initial data, there is a unique weak solution by re- 
ducing the problem to the solution of a Volterra integral equation involving 
only the time variable. It is also shown that the evolution operator uniquely 
extends to a unitary operator in L^{R^). 



Dedicated to Sergio Albeverio 



1. Introduction 



We consider the Schrodinger equation in with an interaction supported 
by n points which are moving on preassigned smooth paths. More precisely let 
a — (ai, . . . be a vector in i?" and let y{t) = {yi{t), . . . ,yn{t)) be n given 

smooth non intersecting trajectories in i?"^. For t G R, let H^ yf^t^f be the Schrodinger 
operator in L^(R^) with point interactions supported at y{t) and with strength a. 
We recall below the explicit definition of Ha^y{t)- 
We are interested in the non-autonomous evolution problem 



(1) ^ = H^Mt)^s{t), Ms)^f 

where s is an arbitrary initial time and / is some (possibly smooth) initial datum. 
An existence theorem for the solution of problem (1) cannot be given using the 
standard theory of non-autonomous evolution problems (see e.g. because of 
the strong dependence on time of the operator domain, and in fact even of the 
form domain of Ha,y(t)- Note that the case of point interactions at fixed positions 
with time-dependent strengths (see |Q , |§ ) is less singular since the form domain is 
constant. 

As we shall see, problem (1) does not have solutions in a strong sense. The reason 
is that, even for very smooth initial datum, the solution exhibits an additional 
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singularity at the position of the moving points and then it does not belong to the 
operator domain. 

In a previous paper ([Q) we studied the corresponding problem for the heat equa- 
tion. For each uq £ D(Ha^y(o)) we proved existence and uniqueness of a strong 
solution, i.e. of a function u(t) belonging for each i > to D{Ha,y(t)), satisfying in 
the L^-sense the equation 



du{t) 

(2) = -Ha^y!^t)U{t), u(0) = Uq 

The proof exploited the smoothing properties of the heat kernel and it cannot be 
generalized to the Schrodinger case. 

In this paper we show that, when interpreted in a suitable weak sense, problem (1) 
has a unique solution. 

More precisely, let Byi^t){'i ') be the bilinear form associated to H^^yi^t) and let Vt be 
its domain (which depends on y{t)). Let 



(3) C^^,)^C^(R^\{y{t)}) 
and notice that C^^^ C D{Ha y^i)) (see (5)). 

We shall prove that for all / G ^'^(s) there is a unique solution of the equation 

(4) ^ («(i), = By^,Mt)Mt)), U^) - / 

for all v{t) e Vt- Moreover V's(i) has a natural representation (see (14)). 
The maps f ^ ips{t), s,t £ R, are isometrics and extend by continuity to unitary 
maps U{t, s) in L^{R^). The maps U{t, s) are continuous in s, t in the strong oper- 
ator topology, and therefore define a time-dependent dynamical system in L^{R^), 
with generator Ha,y{t} at time t. 

Notice that due to the assumptions on the initial data we do not define a flow on 
Vt- We conjecture however that indeed problem (4) defines a flow in Vt, continuous 
with respect to the Banach topology defined on Vt by the bihnear form By(^t)- 
We consider the solution of problem (4) as the first step in the study of the motion 
of a quantum particle (e.g. a neutron) in a fluid, regarded as an assembly of n 
classical particles, each of which acts through a potential of very short range and 
therefore can be considered as a point interaction. The limit n going to infinity for 
the case of the heat equation was studied by us in Q . 

The results presented here are also a preliminary step in the analysis of a class of 
nonlinear models in which the motion of the n classical particles is not preassigned 
but rather determined by the interaction with the quantum particle. 



2. Definitions, motivations and statement of the results 



We have denoted by Ha^y(t) the Schrodinger operator in L'^{R^) with point 
interactions of strength a = (ai, . . . ,a„) placed on the points with coordinates 
y{t) — {yi{t),... ,yn{t)). For the sake of simphcity we have assumed that the 
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strengths are constant, and we shall omit them in the labels from now on. The 
extension to the case where also the strength of the interactions depends on time 
is straightforward, since this dependence on time does not alter the form domain. 
The operator ify(t) is self-adjoint, bounded below, with domain and action given 
respectively for each value of t by 



= I u{t) e L2(i?3) I u{t) = m + qk{t)G{- - yk{t)), 

[ k=l 

m e A0(i) e L^R"), q,it), . . . ,q^{t) e C, 



(5) lim [u{x,t) ~ qk(t)G{x - yk{t)] = akqk{t), k = 1, . . .n 

(6) H^MtM*) = 

Here H'^{R^) is the standard Sobolev space, C denotes the set of complex numbers 
and G is the Green's function 

(7) G{x - x') = (-A)-i(x - x') = / 

It is clear from (5) that the operator domain consists of functions with a regular 

part plus the "potential" produced by the "point charges" qk{t)- The limit in 

(5) is regarded as a boundary condition satisfied by u{t) at y{t). 

We refer to for a complete analysis of this kind of hamiltonians. 

Denote by J-y(^t)jD{!Fy(^t)) the closed and bounded below quadratic form associated 

to -ffy(t) and let By^t^ the corresponding bilinear form. 

One has (see Q for details) 



D{Ty(t)) = {u{t) e L\R^) I u{t) = <j>{t) + qk{t)G{- - yk{t)), 

fc=i 

(8) cbit)eHUR'), \ym^L\R')} 

p n 

^y{t){u{t))^ / dx|V0(x,t)|2 + Vafc|(jfc(i)|2 

k=i 

n 

(9) - Y G{yk{t) - y,{t))qk{t)q,{t) 
where z denotes the complex conjugate of z e C. 

Notice that iVCj ^ L^{R^)^ and therefore the decomposition in (8) is unique. We 
also emphasize that D{!Fy(^t-^) is strictly larger than the form domain of the laplacian. 
To simplify the notation we denote by Vt the Hilbert space D{J-y(t)) equipped with 
the scalar product 

(10) <v{t),u{t) >= By^tMt),u{t)) + f3{v{t),u{t)) 
where f3 > — inf (T(iJj,(()). 

We also introduce the dual space V^* of Vt with respect to the L^-scalar product 
and denote by {^{t),ri{t)) the corresponding duality, ^{t) G Vt, ri{t) € V^* . 
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Finally we define the following set of smooth curves in 



M = {y= {yi,... ,yn) Ivj : R , yj is of class C 
(11) i = 1, . . . n, inf inf \yj {t) - yi{t)\ > a > 0} 

With these notation our main results are the following. 



Theorem 2.1. Let y € M, s € R and f e C^^y Then there exists a unique 
V's(t) €Vt,t€R, such that e V^* and 

(12) *("(^)'^^) =^yitMt),Mt)) Mt)&v, 

(13) Ms)=f 

Moreover ips {t) has the following representation for t > s 

(14) Mt) = Uoit -s)f + iY, dr Uoit - r; • - »(r))g,(T) 
where UQ{t) is the free unitary group defined by the kernel 

(15) Uo{t;x-x') = e'^'{x-x') 



(47rit)3/2 

and the charges qj{t) satisfy the Volterra integral equation 



(16) 
where 



1 /■* 1 f . , , . dB, 



(17) + •'•'^ 



2(a - r) 

^''^ ^^"'^''"^ = ^(^^v^ io ' 

^^^^ ^^■'(^-) = Mbi 



dz e' 
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[20) Wji[t,T) = — 2sft^ — ' 

(21) r)ji(t,T) = / da^=U^{a - T;yj{a) - yi{T)) 

TI" Jt \/t.- O 

A similar representation for the solution holds for t < s (see Sect. 4)- 

Using the representation of the solution we can moreover prove 

Theorem 2.2. The map f '4's{t), s,t € R, extends uniquely to a unitary 
map U{t,s) in L'^{E?). 

The conditions we impose on the smoothness of the curves are not optimal. Optimal 
conditions can be found analyzing in detail the representation of the solution. We 
do not discuss further this problem here but notice that it may be relevant in the 
coupled case. 



3. Some auxiliary lemmas 



We shall construct the solution of (12), (13) for t> s. The case i < s is obtained 
following the same steps and it is outlined in Sect. 4. We start considering Vs(i) 
given by (14) for some functions qj{t). 

In the following we shall drop the dependence on the initial time s. 
We show first that if qj{t) and y{t) are sufficiently smooth, e.g. y G M and 
qj £ Wi^^ (R), then ^p{t) belongs to the form domain Vt- It will also be clear that 
'ijj{t) does not belong to the operator domain even for an arbitrarly smooth qj{t). 
In the second step, using the representation (14) for tp{t), we reduce the solution 
of (12), (13) to an integro-diffcrcntial equation for qj{t). 

In the third step wc show that the resulting equation is in fact equivalent to the 
integral equation (16), which has a unique solution with the required regularity. 
The first result is summarized in the following lemma 

Lemma 3.1. Assume y e M and q^ € Wl^^{R), with qj{s) = and f e C^^y 
Then ip{t) G Vt, where tp{t) is given by (14)- 

Proof. Expression (14) has a simpler form in the Fourier space 

n „t 

(22) ^p{k,t) = e-^'=' (*-«)/>) + J2 J rfre-''='(*--)e-''=-^^Wg,-(r) 
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We prove first that tl>{t) e L^{R^). Due to the regularity assumptions on y, qj and 
/, it is sufficient to prove that 



(23) 
where 



/ dk\Q{k, 

J\k\>l 



t)\^ < OO 



(24) Q{k, t) = dre-'=^(*--)e-^'=-''^(^)9,(r) 

An integration by parts yields 



(25) 



+ 



(27r)3/2fc2 



drk ■ yj(r)e-''='(*-^'-''=-^-'(^)gj(r) + 



-ik-yj{t) 



(27r)3/2fc2 



The only delicate term in r.h.s. of (25) is the second. The explicit computation of 
its Z/^-norm gives 



1 



1 



(2^)3 Jfi3 k^ 



t i-t 



dk—r I dr 



f daqj{T)qj{a)k-yj{T)k-yj{a) 

J s 



xe 



-ik^{(7-T)-ik-{yj{T)-yj{a)) 



dr / daqj{T)qj{a)\yj{T)\\yj{a)\ / ^17(6*, </>) cos^J cos^J 



POO 

(26) X / f^fce"*'^^^'^"'^^"*'^!''^'^'^^"^^'^'^^! 
Jo 

where we have denoted by the angle between k and yj{v). The last integral in 
(26) can be written as 



f 

Jo 



(^j^g-»fe^(o---r)-ife|i/3(r)-j/j(<7)|cose 



(27) 



a -T 



2^/^ 



and similarly for r > cr. 

Using (27) one easily sees that the l.h.s. of (26) is finite and hence one concludes 

that tpit) G L'^{R^). 

Now we have to show that 



(28) |V0(i)| ^ 

Prom (22), (25) we have 
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1 " /■* 



t) — e 



-ik'^ (t — T) — ik-yj (t) 



(27r)3/2fc2 ^ 

\ ' 3 = 1 

(29) =0i(fc,i) + ^2(fc,i) + 03(fc,t) 

The smoothness of / guarantees that |V(/)i(t)| G L^{R^). Concerning (t>2{t) we have 



dkk^\4>2{k,t)\'^ 

< csup / dkk^ 
ft rt 



1 /■* 

^ drg, (r)e 



ik^ {t—T) — ik-yj{T) 



1 



(30) =csnpj dT J daqj{T)qj{cr) j dfc— e^^'^' 

The last integral can be explicitely computed. Using spherical coordinates and the 
position (20), for a > t one has 



27r 



I I . I dp sin {2wji{a,T)p) 

Vcr - TWjl(<7,T) Jq P 



(31) 



2^3/2 



\pL\Ja — T ^ 

(See e.g. Q). An analogous computation holds for a <t. The function Bji{t,s), 
t > has been defined in (19) and it is continuous in both variables and differen- 
tiable in the second one. 

From (30) and (31) one easily gets the estimate for ip2- 

It remains to estimate (pzit)- A further integration by parts yields 



. n t 

{2nf/^,{k,t) ^ -Y^ / drg,(r)fc.y,(r)e-'=^(*-^)-*^^M 
j=i 

n „t 

+ pE / drg,(r)fc.j/,(r)e-'=^(*-)-'=-^^(^) 
+ fc4E / rfrg,(r)(A:.y,(r))^-'=^(*-)-^-«^M 



(32) 



The only delicate term in the r.h.s. of (32) is the third one. Proceeding as in (30), 
one easily sees that its gradient has a finite L^-norm and this concludes the proof 
of the lemma. 
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Remeirk. From (29), (32) we get the following representation for ^{k,t) 



^ ' .5 = 1 



Assuming further regularity on y. and using again integration by parts one easily 
sees that Ax G L'^{B?) which implies Ac/) ^ L'^{E?). 

This means that ip{t) given by (22) does not belong to D{Hy(^t)), i.e. problem (1) 
does not have strong solutions. 

In the next lemma we reduce the evolution problem to the solution of an integro- 
differential equation for qj{t). 



Lemma 3.2. Assume y e M, qj G Wl^^{R), with qj{s) = 0, and f e C^^y 
Then il){t) given by (14) solves problem (12) (13) if qj{t) solves the equation 



4n{Uo{t-s)f){y,it))='l7ra,q,{t)- ^^^^ 



,JTr^jMt)-m{t)\ 

1 /■* , . . ^Bii(t,T) Vi ^ f* , , ^^^(<,r) 

Vi^jr^Js y/t-T V^^'frtJs y/t-T 

(33) 

where Bji{t,T) and Aji{t,T) are given in (19), (18). 

Proof. Prom lemma 3.1 we know that i/'(i) G Vt and then the r.h.s. of (12) 
is well defined. Now we check that ^^^^ G V^*. A direct computation yields (see 
(22),(25)) 



- (^^^e/ («iW^-"-^^0 



= -ik'^(j){k,t) 



(34) 

For v{t) G Vt, we write 



(35) y^t)=C{t) + Y^q]{t)G{--y,{t)), C{t)GH\R') 



Using (34) we have 
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(36) (rW, ^) = -ij^^ dxVe{x,t) ■ Vcl>{x,t) 

which is obviously finite. Moreover 



dk 



dk- 



dt J Jji3 (27r)3/2A;2 dt 

(27r)3/2fc2 



fl3 



^ ft J 



= -!{l'o(t -»)/)(»(*)) 

" l-t 



-ik^{t-T)+ik-(vj(t)-Vl(T)) 



{t)-Vl{T)) 



The last integral in the fc-variable can be explicitely computed. We introduce 
spherical coordinates k = {r,9,(j)) with polar axis directed along yj{t) — yiij) and 
yi{'^) = (|y(('^)l)^)0). Using the formula 



(38) 

we have 



k-yiir) =r\yi{T)\ (■ 



cos cos ^ + sin ^ sin 6 cos < 



L 



kj_yi{T)_ ik^(^t-T)+ik-iy,it)-yii 

47r|y;(T)| cos^ 



Wt-T\yj{t)-yi{T) 



dpe cos{2w j I (t,T)p) 



1 



oo -ip 



dp sin{2w ji{t, t)p) 



(39) 



2wji(t,T) Jo p 

2(^)3/2|y;(r)|cos^ 
iVi^/t-T\yj(t) - yi{T)\ 

4(7r)3/2|y;(T)|cos^ 
ViVt-T\yj{t) - yi{T)\ Wji{t 



1 



Wjl{t,T) Jo 



dve"' 



{t,T)l 



Wjl{t,T) 



dvv e 



^^i^/T■ 



-.Aji{t,T) 



where the function Aji{t,T) has been defined in (18). From (37), (39) we find 



10 



G.F. DELL'ANTONIO, R. FIGARI, AND A. TETA 



- yjit)), ^) = ~^{Uo{t - s)f){y,{t) 

(40) 

which implies e V^^*. Using (36), (40), it is now easy to check that the evolution 
equation (12) reduces to the integro-difFerential equation (33) which is satisfied by 
hypotheses. This concludes the proof of lemma 3.2. 



4. Proof of theorem 2.1 




We shall use the results of lemma 3.1 and 3.2 to complete the proof of theorem 2.1. 
Let us fix y S and the initial datum / e ^'^{s)- -^^"^ t > s we consider the 
equation (16) for qj[t). 

It is a Volterra integral equation containing the Abel operator 



(41) {Lq,) (t) 

and the integral operators 



(42) {C,q,){t)^ I dTq,{T)C,{t,T) 

t 

(43) {Djiqi)= I dTqi{T)D,i{t,T) J^l 
The datum of the equation 



(44) h,{t)^^ / dT 



f\JUo{T-s)f){y,{T)) 



is the result of the application of the Abel operator L to the function 

47r(C/o(i ~ ^)f){yj{t))- Due to the smoothness of / it is obviously true that hj e 

Wl^^iR) and/i,(s) = 0. 

By direct inspection of (17)-(20), one verifies that ii y E A4 then Cj{t, r) is contin- 
uous in both variables and it is differentiable as a function of t. 
By a detailed analysis of the expression (21) (which we omit for brevity), one can 
also show that Dji is a bounded operator in Wi'{R). 

The same is true for the Abel operator (see e.g. [g) and we conclude that equation 
(16) has a unique solution qj G W^^^^R), with qj{s) — 0. 

Now we apply the Abel operator to equation (16), and make use of the fact that 
for rj differentiable with r]{s) — one has 
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(45) 1[(L)2^]=,^, ±[Lr,){t) = {Lm 

The resulting equation reads 



d d 
-i{Lqj){t) + 47ra,«,(t) - i-{LCjq){t) - i J2 (^^m) (0 

(46) =4n{Uo{t-s)f){yj{t)) 

The integral operator Cj can be rewritten as 

= -i f da^ r drq.ir) ('-^^ + ^ f 



(47) +1 r*,^ r*,,M2«i5ii^ 

Using again the first equation in (45) we have 

(48) / (ir- 



Concerning the integral operator Dji we have 

(LDjiqi) (t) = drqiir) J daUo{a - r; yj{a) - yi{T)) 

• /' /" ,r,M?pd ^^f^r dr,,ir)^ 
(49) -i [ daqi{a) / 

If we substitute (49), (48) into equation (46) we find that the charges qj{t) satisfy 
the integro-differential equation (33). By lemma 3.2 this means that if the qj{t) 
solve equation (16) then ip{t) given in (14) solves the evolution problem (12), (13) 
for t > s. 
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We now briefly consider the case of the backward evolution. More precisely, given 
the initial time i e i? and g G C^t)i we want to find € K,, for s < t, satisfying 

the equation 

(50) Mt) = g 

Again we start representing iptis) as 

(51) Ms) - t/o (i ^s)9-^Y. "^^^oir -s;-- y,{r))q,{T) 

1=1 

for some functions qj and then we determine qj in such a way that (51) solves (50). 
The steps are similar to the case of the forward evolution and will be omitted. We 
only write the integral equation which is satisfied by qj 



q, is) + a, ^ /* + / drq, {t)C, (r, s) 



'I 



1=1,1^3 

(52) 



[ dTqi{r)D,i{T,s)^^ f dr 



f\jU^it-r)9){yj{T)) 



Finally the uniqueness of the solution of problem (12), (13) easily follows from the 
fact that for any solution of (12) the L^-norm is conserved. 



5. Unitary evolution 



In this section we give the proof of theorem 2.2 following the idea developped 
in for the case of point interactions at fixed positions with time-dependent 
strengths. 

We fix s, f e i? and, without loss of generality, we take s < t. 

By theorem 2.1 we have existence and uniqueness of the forward evolution ips{t) 
and the backward evolution iptis) for smooth initial data, denoted respectively by 
/ and g. 

Moreover the linear maps 



(53) f^Mt), 9^ Ms) 

are both defined on a dense set of L^{R^) and are isometrics. Then they can be 

uniquely extended to isometrics on L?{B?). 

We shall denote them respectively by U{t,s) and U{s,t). 

In order to prove that they are unitary maps, we have to show that the adjoints 
U*{t,s) and U*{s,t) are also isometrics. This fact will follow from the equalities 



SCHRODINGER EQUATION WITH MOVING POINT INTERACTIONS IN THREE DIMENSIONS 



(54) U*{t,s) = U{s,t), U*{s,t) = U{t,s) 

Here we prove the first equality in (54) (the second is obtained in the same way). 

For the sake; of clarity, it is convenient to rewrite U {t. s) using various integral 
operators. In particular we use the Abel operator L on L'^{[s,t\) defined in (41) 
and its adjoint L*, the operator T : L'^{R^) ®'^^^L'^{[s,t]) defined by 

(55) {Th)j{T) = / dxUoir; x - yj{T))h{x) 
and its adjoint T* . Moreover, from (17), (21), we have 

(56) Cj=LRj, Dji=LSji 
where Rj and Sji are the integral operators defined by 

^""^'^^^^^ =-7rJs '"7^ r^-^-^"' '^^ + ^ 2(V-.) ) 

(57) 



(58) {Sjiqj){T) = daqj{a)Uo{T~cr;yj{T)-yi{a)) 

We also introduce the corresponding operators on the space of vector valued func- 
tions q(t) = (gi(t), . . . , g„(i)), i.e. 

(59) {Rci)j{r) = iRjqj)ir) 

n 

(60) (^q)i(T)= iSMi)ir) 

(61) {Aci)j{t) = 4TTa,q,iT) 

Using the above notation we rewrite equation (16) in the form 

(62) (i + L{A + R + S)q = 2y^LTUS{s)f 

From (14), (55) and (62) we obtain the following representation of U{t,s)f, for 
s < i and / e C^,) 

(63) U{t, s)f = Uo{t - s)f + iUo{t)T* [I + L{A + R + S)]'^ LTU^{s)f 

The same procedure can also be applied to obtain the representation of U{s, t), for 
s<tandge C^^^ 

(64) U{s, t)g = U^{t - s)g - iUo{s)T* [I + L*{A + R* + S*)]'^ L*TU^{t)g 
A straightforward computation gives 
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iU*{t,s)gJ) = {g,U{t,s)f) 

= {U;{t - s)g, f) + i (Uois)T*L* [{I + L{A + R + S))-'] * TU^{t)g, /) 
= {U^{t - s)g, /) + I {Uois)T* [I + L*{A + R* + S*)]-' L*TUS{t)g, /) 

(65) 

Thus, from (64), (65), we have U*{t, s) = U{s, t) on a dense set and then on L'^{B?) 
and this concludes the proof of the theorem. 



Acknowledgements. One of us (G.F.D.A.) is grateful to the Von Humboldt 
Foundation for the generous support. 



References 

[1] Albeverio S., Gesztesy F., Hogh-Krohn R. and Holden H., Solvable Models in Quantum 

Mechanics, Springe- Vcrlag, New- York, 1988. 
[2] DcU'Antonio G.F., Figari R. and Tela A., A limit evolution problem for time-dependent point 

interactions, J. Func. Anal, 142, 1996, 249-275. 
[3] Deir Antonio G.F., Figari R. and Teta A., Diffusion of a particle in presence of N moving 

point sources, Ann. Inst. H. Poincare A, 69, 1998, 413-424. 
[4] Erdely et al.. Tables of Integral Transforms, vol. 1, McGraw-Hill, New- York, 1954. 
[5] Gorenflo R. and Vessella S., Abel Integral Equation, Springer- Verlag, Berlin-Heidelberg, 1991. 
[6] Sayapova M.R. and Yafaev D.R., The evolution operator for time-dependent potentials of 

zero radius, Proc. Stek. Inst. Math., 2, 1984, 173-180. 
[7] Tcta A, Quadratic forms for singular perturbation of the laplacian, Publ. R.I. M.S., 26, 1990, 

803-817. 

[8] Yafaev D.R., Scattering theory for time-dependent zero-range potentials, Ann. Inst. H. 

Poincare' A, 40, 1984, 343-359. 
[9] Yajima K., Existence of solutions for Schrodinger evolution 

equations. Comm. Math. Phys., 110, 1987, 415-426. 

DiPARTIMENTO DI MATEMATICA, UnIVERSITA' DI ROMA "LA SAPIENZA" , ITALY AND LABORA- 
TORIO INTERDISCIPLINARE, S.I.S.S.A., TRIESTE, ITALY. 

E-mail address: dellantonioSmat.uiiiromal.it 

DiPARTIMENTO DI SCIENZE FiSICHE, UNIVERSITA' DI NAPOLI, ItALY. 

E-mail address: figariSna. infn.it 



DiPARTIMENTO DI MATEMATICA, UNIVERSITA' DI ROMA "LA SAPIENZA" , ITALY. 

E-mail address: teta9mat.uniromal.it 



